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■ If one assumes a translationally invariant motion of the nucleons relative to the c. m. 

position in single particle mean fields a correlated single particle picture of the nuclear wave 
function emerges. A single particle product ansatz leads for that Hamiltonian to nonlinear 
equations for the single particle wave functions. In contrast to a standard not translationally 
invariant shell model picture those single particle s-, p- etc states are coupled. The strength 
CO ' of the resulting coupling is an open question. The Schrodinger equation for that Hamiltonian 

can be solved by few- and many -body techniques, which will allow to check the validity or 
non-validity of a single particle product ansatz. 

Realistic nuclear wave functions exhibit repulsive 2-body short range correlations. There- 
fore a translationally invariant single particle picture - if useful at all - can only be expected 
beyond those ranges. Since exact A = 3 and 4 nucleon ground state wave functions and be- 
' Qjj ■ yond based on modern nuclear forces are available, the translationally invariant shell model 

^ ' picture can be optimized by an adjustment to the exact wave function and its validity or 

<**] \ non- validity decided. 

| §1. Introduction 

The shell model for the nucleus has a long tradition. However, in its standard 
form expressed in single particle variables it is plagued by violating translational 
fSj ■ invariance. Various methods have been suggested to remedy this situation, like for 

instance the generator coordinate methodP'^'®'® Clearly, if the shell model is 
realistic at all, the motion of the individual nucleons in a mean field happens in a 
translationally invariant manner, namely as a function of Ui = Xj — X, where X{ are 
the individual coordinates of particle i and X is the c. m. coordinate. However, this 



X 



set of coordinate vectors u« obeys the obvious condition Ylt=i ™i = 0> correlating the 
motion of all particles. 

It is the aim of the present investigation to work out the consequences of choosing 
the coordinates U{ for a shell model picture. 

In Section II we provide some formal basis for this specific choice of coordinates. 
We restrict ourselves in this first investigation to systems of three and four nucleons. 
Furthermore, our most simplistic ansatz for the wave function shifts the antisym- 
metry requirement to the spin-isospin space, which leads to a symmetric space part 
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under particle permutations. Then the very first ansatz for the space part is 

A 

<p(ui)=HR( Ui ), (i-i) 

i=i 

with A = 3 or 4. Here only s-wave states are assumed. We denote such a form a 
correlated single particle picture. 

The nonlinear equations for 3 and 4 particles for the state R{u) assuming a sum 
of single particle potentials, V = Ylt=x V( u i)i are presented in Section III. 

In the case of the harmonic oscillator potential the nonlinear equations can be 
solved analytically and it is shown that the ansatz is indeed the correct one. 

Obviously the question arises whether the ansatz (|l-ip for the wave function 
is at all valid for general single particle mean field potentials. To that effect the 
Hamiltonian, including the sum of single particle potentials, can be expressed in 
standard Jacobi variables. This is displayed in Section IV. 

The resulting Schrodinger equations for 3 and 4 particles (in this case bosons 
for the space part) can be solved exactly in the form of the Faddeev-Yakubovsky 
equations, which will be formulated in that Section. Having the exact wave function 
at ones disposal, one can then investigate how well the above shell model ansatz 
is realized, whether contributions beyond s-wave are needed, or whether that hope 
is not realistic at all. An optimization algorithm relating the shell model ansatz to 
the exact wave function is presented in Section V. 

The numerical investigations for solving the nonlinear equations for R(u), for 
solving the Faddeev-Yakubovsky equations for the shell-model Hamiltonian and for 
the optimal extraction of R(u) from the exact wave functions is left to forthcoming 
investigations. 

The main task however remains. The realistic nuclear wave function is deter- 
mined by two- and three- nucleon forces. First estimates in the case of the a-particle 
indicate that even small contributions from proper 4N forces^® are needed. Based 
on these forces numerically exact wave functions are nowadays routinely generated 
for three and four nucleonsP'®'^'^'^'^'^ The question however arises, how 
well these wave functions for pair distances larger than a certain distance ro can be 
represented in the form of a correlated shell model ansatz like in (|l-ip . or whether 
higher partial waves and more complicated symmetries with respect to space-, spin- 
and isospin parts of the wave function are required. Clearly for pair distances smaller 
than ro short range repulsive features are present in the realistic wave functions which 
can not be represented in the shell model form. On the other hand it is known that 
the short pair distance behavior is essentially universal for light nucleP^'^'^ aside 
from proper normalization, which might allow an overall description: short range 
depletion and correlated shell model feature at larger distances. The value ro is 
expected to be somewhat smaller than 1 fm. Section VI provides some suggestion 
on how an optimal extraction of a correlated single particle picture can be obtained 
from realistic three- and four-nucleon wave functions. We summarize in Section VII. 
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§2. The Formal Basis 

The translationally invariant single particle coordinates for n particles are de- 
fined as 



n n 

- -> n ~ - - 

Ui = Xi / Xj — X{ / Xj [Z'li 

n. ± — ' J n n £ — ' J 



n < — ' n n 



Since Y^i=i = 0> the mapping from the n Xi to the n can not be inverted and 
we choose the new variables as the first (n-1) Hj together with the cm. coordinate 
X 



1 n 



n 
3=1 

It is a straightforward exercise to express the kinetic energy in terms of the new 
variables 

2m fc 2m n n 2mn 

k=l i=l i^j 

Clearly, the first part refers to the relative motion, the second part to the c. m. 
motion. While the choice of Jacobi coordinates avoids mixed terms in the kinetic 
energy, here they are unavoidable. 

Let us now restrict ourselves to three and four particles. If one chooses a Slater 
determinant with equal space-dependent single particle wave functions, the symmet- 
ric part of the form (|1T|) factors out and one is left with a totally antisymmetric 
spin-isospin part. For a proton-proton- neutron (ppn) system this has the form 

X3 = l(* = o|)r = !>!(« = i|)fir = I > 

-\(t = l\)T= l ->\( S = U l -)S= l ->, (2-4) 

where the two-body spin or isospin state is coupled with the spin or isospin | of the 
third particle to total spin S = \ or total isospin T = h. This state together with a 
symmetric space part is knowrP^'^ as the principal S-state for realistic 3 He wave 
functions and contributes with about 90% to the norm. This result by itself clearly 
indicates that this choice of the one Slater determinant can not exhaust the full wave 
function but at least a very large portion of it. 

For the ppnn system the totally antisymmetric spin-isospin part of the wave 
function has the form 

X4 = (1 - P23 - ft4)[(~)0(~)05 = > 

(1 + ^13^24) |(~)H >12 |(~)1 " 1 >34, (2-5) 
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where Pjj is a transposition of particles i and j. That state has total spin 5 = 
and total isospin T = 0. Again in relation to the norm of a realistic a-particle wave 
function it accounts for about 90% PJ 

Now we provide some formal properties, whose verification is left to the reader. 
The Heisenberg commutation relations 

[Wka,Ujp] = hj-T&aP (2-6) 

are obeyed, where Wk a = gjjp are components of the conjugate momenta. 
The relative orbital angular momentum has the form 

L re i = ui x tV Mi + u 2 x -V U2 , (2-7) 

i i 

which justifies that standard Clebsch-Gordon coupling in the variables u\ and u 2 
can be used. 

Using (|2-3|) for n = 3, the translationally invariant shell model Hamiltonian is 
given by 

H 3 = + V 2 U2 - V U1 • V tt2 ) + + V(u 2 ) + V(u 3 ), (2-8) 

where U3 = |«i + v,2\. Obviously, a separation of variables is not possible. However, 
the symmetry of the kinetic energy under particle exchanges is valid: 

+ v 2 U2 - V W1 • v U2 = v 2 U2 + v 2 U3 - V U2 ■ V U3 

= V^ 3 + V^-V U3 -V U1 . (2-9) 

In the case of four particles the translationally invariant shell model Hamiltonian 
is given as 

H4 = -^(v^ + v^ 2 + V 2 U , A - |(V U1 • V U2 + V M1 • V U3 + V U2 • V U3 )) 

+ V( Ul ) + V(u 2 ) + V(u 3 ) + V(u 4 ), (2-10) 

with U4 = I iTi + U2 + u,3\- All the formal relations corresponding to (|2-6p .( [2^7]) . and 
(|2-9p are valid for four particles as is expected. 

§3. Nonlinear Equations for the Translationally Invariant Shell Model 

States 



For three particles the simplest ansatz for a symmetrical space part wave function 

is 

$(m,u 2 , u 3 ) = R{ui)R(u 2 )R(u 3 ), (3-1) 

with U3 = \u\ + u 2 \. 

It is straightforward, though tedious, to evaluate the action of the kinetic energy 
in (|2-8[) onto (|3-1|) . If we put R(u) = the Schrodinger equation based on #3 
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and the ansatz (|3-ip results in 
1 



where 



L r (u\)r{u 2 )r{u 3 ) + r{u\)r (u 2 )r(u 3 ) + r(ui)r{u 2 )r (u 3 ) 
6m 

- (r (ui) )r(u 2 )(r (u 3 ) )u\ ■ u 3 

U\ u 3 

i \! '( \ r ( n 2)w s r(u 3 ) A 

- r(ui)(r (u 2 ) )(r (u 3 ) )u 2 • ii 3 

- (r (ui) )(r (it 2 ) )r(U3)«l ' H 

ni u 2 

+ (V(ui) + F(u 2 ) + F(u 3 ) - E)r{ Ul )r{u 2 r{u 3 ) = 0, (3-2) 



U\ + Ul- u 2 

Ul -u 3 = -■ 



U3 

„ A M 2 + U 2 ' Ul . . 

"2 • ^3 = (3-3) 

Here the independent variables are u±,u 2 and x = u\ ■ u 2 . 

We can not exclude that higher partial waves should be included. The simplest 
ansatz for a p-wave admixture is given by 

$i(ui,u 2 ,u 3 ) 

= Ri(ui)Ri(u 2 )R(u 3 )ui ■ u 2 + R(ui)Ri(u 2 )Ri(u 3 )u 2 ■ u 3 

+ i?i(ui)i?(u 2 )i?i(u 3 )ui • u 3 , (3-4) 

where for the sake of simplicity we assumed that the 'third' state, which is not 
involved in the p-wave admixture, remains unchanged. We leave it to the reader to 
derive the resulting equation. 

For four particles the most simple ansatz is 

$(ui, u 2 , u 3 , u 4 ) = R(ui)R(u 2 )R(u 3 )R(u^), (3-5) 

with \u4 = \ui + u 2 + u 3 \. The resulting equation based on H$ is 
3 // 

- — [r {ux)r(u 2 )r(u 3 )r(u4) 

+ r(ui)r(u 2 )r(u 3 )r (u 4 ) + r(u\)r (u 2 )r(u 3 )r(u4) 
+ r(ui)r(u 2 )r(u 3 )r (u 4 ) + r(ui)r(u 2 )r (u 3 )r(ui) 
+ r(ui)r(u 2 )r(u 3 )r (u 4 ) 

-2(/(u 1 )-^)r(t l2 )rW(r'W 

-2r(u 1 )(r'(u 2 )- r ^-)r(u 3 )(r , (i M ) 
u 2 

- 2r(ui)r{u 2 )(r \u 3 ) )(r'(u<l) 

dm 



r(it 4 ). 








■ U 4 


M 4 






r(u 4 ) 








• U 4 


u 4 






r(u 4 ) 






)«3 


■ u 4 


u 4 
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+ (r (m) )(r (U2) )r(u 3 )r{u4)ui ■ u 2 

Ui U 2 

- ((r( Ul )(r'(u 2 ) - r -M)r(u 3 )(r'(u 4 ) - ^) 

U 2 U4 

+ r( Ul )r( U2 )(r>(u 3 ) - ^)(r'(u 4 ) - ^) 

U 3 U4 

+ r( Ul )(r'(u 2 ) - ^-) r (u 3 )(r'(u 4 ) - ^)))u 2 • u 4 

- (((r'(m) - r -M) r ( U2 )r(u 3 )(r'(u 4 ) - ^) 

+ (r'(«i) - ^)r(u 2 )r(u 3 )(r'(u 4 ) - • u 4 

U\ U 4 

. .. . r(ui) . K^ 3 ) \ / \ ~ 
+ (r (tti) )r(u 2 )(r [u 3 ) )r(u 4 )«i • u 3 

- r(«i)r(« 2 )(r (it 3 ) )(r (w 4 ) )u 3 • u 4 

ti 3 u 4 

+ r(ui)(r (u 2 ) )(r {u 3 ) )r(ii 4 )u 2 • u 3 

u 2 u 3 

+ (V(i*i) + ^(7/2) + V(« 3 ) + ^(u 4 ))r(ni)r(n 2 )r(n 3 )r(ti 4 ) 

= £ , r(wi)r(u2)r(u 3 )r(n 4 ). (3'6) 

Again, extensions to the ansatz ()3-5[) are obvious. Already the presence of the 
explicit angular dependence in (|3-2|) and (|3-6p suggest that higher orbital angular 
momentum admixtures are likely and that the most simple ansatz for the ground 
state may be poor. 

Choosing the mean field potential V(ui) to be a harmonic oscillator, V(ui) = 
^jrU 2 , the nonlinear equations, (|3-2p and (|3-6p . can be solved analytically. As exam- 
ple we consider four particles and introduce standard Jacobi coordinates 

X = x 2 - x 3 

y = xi- -{x 2 + x 3 ) 
„ 1,-. 

Z = X4- -{Xi + x 2 +x 3 ) 

X = -{x\ + x 2 + x 3 + f 4 ). (3-7) 
Then the potential energy V = ^(^x 2 + |y 2 + \z 2 ) as well as the kinetic energy 

r„, = -^(2v= + |vj + ivl) (3-8) 

allow for a separation of the variables with the result 

${x, y, z) = e -^ x2 +b 2 +h 2 ) = Y[ R( Ui ). (3-9) 

where 

R{u) = e~— u . (3-10) 
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This goes with the lowest energy E = |u>. The corresponding result for three 
particles, now for E = 3w, is 

3 

${x,y) = Y[R(ui), (3-11) 

i=i 

with the same function R(u). It is straightforward to verify that (|3-9p and (|3T1|) 
fulfill the nonlinear equations (|3-6p and (|3-2p . 

§4. The Faddeev-Yakubovsky Equations for Three and Four Particles 

The shell model Hamiltonians H3 and H4, Eqs. (|2-8p and (|2-10p . can be rewritten 
in terms of standard Jacobi coordinates. This allows one to solve the two Schrodinger 
equations exactly in the form of the Faddeev-Yakubovsky equations and therefore 
to test the quality of shell model ansatz. For three particles one defines the Jacobi 
coordinates as 

X = U2 -u 3 = 2u 2 + Ul 
1 3 

y = ui - -{u 2 + U3) = 2«i> I 4 ' 1 ) 

or 

Ul = -y 

1 -. 1 -. / \ 

This gives for the Hamiltonian 

H * - ~i v i - i v » + + v «¥ - 5*» + v "¥ + 5*» (4 ' 3) 

The above expression has a formal similarity to a three-body Hamiltonian composed 
of two-body forces: 

#3,26 = ~ + V 2b (x) + V(|~a* + y\) + - y\) . (4-4) 

However, Eqs. (|4-2p and (|4-4p are different. Nevertheless the formal structure of the 
Faddeev equation^ can be used. The three-body bound state obeys 

3 3 

# = 6?oXVi* = Z>«' (4-5) 

i=l i=l 

where Go represents the free three-body propagator, and Vi = V(ui). Then one 
arrives in a standard manner® at 

A = G Ti^2^, (4-6) 
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where Tj obeys the Lippmann Schwinger equation 

Ti = Vi + ViG Ti (4-7) 

Because of the identity of the particles one arrives at the well known form for the 
total state 

^ = (l + P)Vl, (4-8) 

with P = P12-P23 + P13P23, which is a sum of a cyclical and an anticyclical permu- 
tation of three particles. One Faddeev equation is sufficient, namely 

ih = GqTxPVi (4-9) 

The Faddeev equation can be solved in configuration space as an integro-differential 
equation or, what we prefer, in momentum space as an integral equation. In the lat- 
ter case one needs the momentum space representation of the shell model potential 
as well as of the Lippmann Schwinger equation in terms of the conjugate momenta 
p x and p y of the Jacobi momenta x and y. 

With standard (unit) normalizations it results in 

1 2 11 

< xy\u\u 2 >= (-) 3 5(«i - -y)S(u 2 - -x + -y) (4-10) 

Furthermore, as consequence of the locality assumption 

< ^ x v! 2 \V{ui)\uiU2 >= S(u 2 - u' 2 )5(ui - u'^Viui) (4T1) 
and using (|4-10p one obtains 

< i?J y \V( Ul )\p x py >= 5(p x - x )-Jxs j ^ye^y-^-W^y). (4-12) 

Due to that structure the T-matrix element in ()4-7p must have the form 

2 

< $Jy\T x \p x py >= 6<jf x - p x )htf y ,p y , z = E-^), (4-13) 

where t\ obeys 

h(tf y ,Py,z) = V x (j?y,py) + / d 3 p'yV(j? y ,Py) * — t x (p' y , p y , z) ■ (4~ 14) 

m 4m^y 

For two-body forces the 5-function in (|4-13|) would have been for the spectator 
momentum p y . We assume that the mean field forces are spin- independent and 
require symmetry in the spatial part. 

TjJgjJ) guch a system has been shown to be easily solvable using directly momen- 
tum vectors and thus avoiding any partial wave decomposition. We follow the same 
approach. Then (|4-9p . using ()4-13j) has the form 

If. p 2 

<M#i >= p| _ _g_ 2 J d 3 p y ,t 1 (py,tf y ,Z = E--?-) 

m 4rn,Py 



On a translationally invariant nuclear single particle picture 9 

J d Z p x nd Z p y n < PxPy'\P\p x "P y " >< P x "P y " IV'l > ■ (4'15) 

The permutation matrix element is well known^* and is given as 

80 2 4 4 2 

< p x Py'\P\p x "P y " > = (g) {S(p y > + + -p x ")S(p y// - -p x - -p x ») 

2 4 4 2 

+ S (Py' ~ %Px ~ ^P X ") 6 (Py" + ^Px + gPa,")) ( 4 ' 16 ) 

Therefore, Eq. f)415[) turns into 

1 8 o /" , 

< PxPyWl > = -2 — (5) / d Px" 

F, — t-Z- — 77^ «5 J 

m Ami y 

3 4 p 2 4 2 

-jP* - g^x" ' Z = £ - ^) < Px" > g& + gPx" 1^1 > 

3 4 p 2 4 2 

+ hiPy, -jPx + « = E — -^-) < p x ",--p x - -p x » \ipi > . (4-17) 

Because of the uniqueness of the solution, any solution of (|4-17|) has the property 
< — PxiPylipi >=< Px-.Py^i >■ This equation is can then be solved by iteration 
using a Lanczos type algorithm.^ 

As follows from (|4- 16|) the total state given by fj4-8f) has the form 

< Px,PyW > = < Px,Py\^l > +(jf < ~^Px ~ ^Py,Px + ^Py\lpl > 

3 13 1 

+ (^f < ~^Px + ^Py, ~Px + ^l^l > • (4T8) 

In the case of four particles we use the Yakubovsky equations For four bosons 
and two-body forces this has been solved rigorously the first time inP21> Now we have 
different potentials depending on the relative coordinates Uj, which require a renewed 
derivation. Starting from 

4 4 

«P = G £v(ui)!P = 2^ (4-19) 

i=i i=i 

one arrives in a standard first step at 

■fa =G T 1 (^ 2 + V'3 + V'4), (4-20) 

where T\ obeys the Lippmann Schwinger equation (|4-14p . (Note however, the mod- 
ified free four-body propagator.) 

In the spirit of the Yakubovsky scheme one regards a three- body subsystem by 
defining 

^l;123 = Gori(^2+^3) (4-21) 

and a remaining component 

^l;M = G Ti^ 4 . (4-22) 
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Then 

fa = ^1;123 + V>1;1,4- (4-23) 

Correspondingly one defines 

V>2;231 = G T 2 (lp3 +1pl) 
^3;312 = G T 3 (^i + 7p 2 ) 
^2;2,4 = G T 2 ^ 4 

V>3;3,4 = G T 3 ^ 4 (4-24) 

where 

V>2 = V>2;231 + ^2;2,4 

V>3 = + ^3;3,4- (4-25) 

Then (j¥ZED and ( yield 

^15123 = G Ti(V>2;231 + V>3;312 + ^2;2,4 + ^3;3,4)- (4-26) 

Due to the identity of the particles one has 

^2;231 + ^3;312 = ^1;123 

^3:3,4 = -P23^2;2,4, (4-27) 

and (|4-26p can be rewritten as 

(1 - G T 1 P)^ 1;123 = G Ti(l + P 2 3)^2;2,4. (4-28) 

The left hand side by itself defines a three-body problem. After inversion one obtains 

V>1;123 = G f(l + P 23 )V>2;2,4, (4-29) 

where T obeys 

f = Ti + TxPGqT. (4-30) 
It remains to consider (|4-22|) . which in analogy to (|4-25p has the form 

V>2;2,4 = G T 2 (V'4;413 + V>4;4,2)- (4-31) 

Using now 

^4;4,2 = ^24^252,4 (4-32) 

we rewrite (|4-3ip as 

(1 - G T 2 P 24 )^2;2,4 = G T 2 V4;413. (4-33) 

Inversion yields 

■02;2,4 = G T^ 4; 4i 3 , (4-34) 
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where T obeys 

f = T 2 + T 2 P 24 G f (4-35) 
Finally permutation symmetry yields 

V"4;423 = -Fb-FWlsm (4-36) 
and one ends up with two coupled equations 

■01;123 = G f(l + P 2 3)^2;2,4 

V>2;2,4 = G fP 23 P U i>l;l23- (4-37) 

The total wave function is now given as 

<P = Ipi + 1p 2 + Tp 3 + tpi = (1 + P)V>1;123 + ^14^23(^1; 123 

+ ^2;2,4) + (1 + ^24 + Pl 2 ) V>2;2,4 ■ (4-38) 

While a corresponding coupled set based on two- body forces has been rigorously 
solved^ 1 in a partial wave representation, it is also possible to directly use momentum 
vectors as has been demonstrated inP^ 

We would propose to follow that second option. We leave it to the reader to 
work out the explicit momentum space representation of (|4-37j) and (|4-38p in terms 
of appropriate Jacobi momentum vectors. 

§5. Shell Model Ansatz versus Exact Wave Function 

The solution of the Faddeev equation (|4-17p yields the full three- dimensional 
three-boson Faddeev component in momentum space. This is the input for the full 
wave function given in (|4T8p . Since we search for the lowest energy state, is a 
scalar and therefore depends only on 3 variables 

< Px,P V \& >->" &(Px,Py,Px -Py) (5T) 

As a consequence, the dependence of the configuration space wave function 
< xy\fJf > will also reduce to a three- variable dependence ^{x, y, x ■ y): 

< xy\W > = — ^ J d 3 Px d 3 p y e i( -P x i;+ Pyy ) ^(p x ,Py,p x ■ p y ) 

= 7^A3 j d 3 p x d 3 p y COS(p x ■ X+Py- y)<F{p x ,Py,p X ■ Py) 

= W(x,y,x-y). (5-2) 

We used the reality property of \P to replace the exponential by the cosine. 
The expectation is now that 

VsM&y) = R{ Ul )R(u 2 )R(u 3 ), (5-3) 

with us = \u\ + u 2 \ being a good approximation to \P(x,y,x ■ y). In the case of the 
harmonic oscillator this is exactly fulfilled. 
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In general one faces the task to minimize \&(x,y,x-y) — R(ui)R(u 2 )R(u3)\ for all 
x,y,x-y 01 u\,u 2 ,ui-u 2 . Explicitly this requirement is \\P(x, y,x-y) — R(^y)R(\^x — 
\y\)R(\\x + \y\)\ or \V(\ti x + 2u 2 \, § Ul , "fr+^ffi " + u 2 \) to be 

minimal. 

Instead of an optimized pointwise adjustment one can try an average adjustment 
minimizing 

3 u± + 2u 2 • u\ 



J du\du 2 dui ■ U2(lP(\ui + 2u 2 \, ^iti, 



-i?(wi)i?(« 2 )i?(|ui+« 2 |)) 2 (5-4) 

in relation to the choice of R(u). For instance, one can expand R(u) into harmonic 
oscillator wave functions <P m (u), where ^jjr is optimally adjusted to the given mean 
field potential V(u). 
Thus 

R(u) = £& n (u)C m , (5-5) 
m 

and the set C m is to be varied minimizing the above integral. Differentiating with 
respect to and putting the result to zero yields a nonlinear relation for the coef- 
ficients C m . This might be solved by an iterative procedure allowing first Co ^ 0. 
Then keeping also C\ ^ in addition one might start with Co from the previous 
step and determine C\. Finally one can iterate the nonlinear equation for Co and 
Ci starting with the values found before; etc. Very likely, however, one has to allow 
in addition for p-wave admixtures as given in (|3-4p and possibly even higher orbital 
angular momentum values. 

The direct solution of the nonlinear equation (j3-2|) poses a severe problem. More- 
over, very likely p-wave and possibly higher order admixtures have to be taken into 
account, which requires an extension of the nonlinear equation (|3-2p as mentioned 
above. Discretization in the ui,u 2 ,x -y - values is necessary and iterative procedures 
appear unavoidable. Thereby each run is of course an eigenvalue problem for the 
energy E. 

In the case of four nucleons the symmetric state of lowest energy is again a scalar 
and thus depends on 5 variables: 

# = &{x,y,z,x -y,x ■ z,y ■ z), (5-6) 

where x, y, z are one choice of standard Jacobi coordinates. The optimal extraction 
of R(u) in 

V S m{x, y, z) = R{ Ul )R(u 2 )R(u 3 )R(u 4 ) (5-7) 

and possibly higher angular momentum admixture follows analogous strategies as 
for three nucleons. 

§6. Realistic Three- and Four-Nucleon Wave Functions 

Based on modern nuclear forces fike^'^'GIS combined with three-nucleon (3N) 
forces of the Tucson-Melbourne type^D or based on the most recent consistent two- 
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and three-nucleon forces generated from chiral effective field theory^ numerically 
exact solutions of the Faddeev - Yakubovsky equations are available. If a correlated 
single particle picture applies at all it can only be valid beyond a certain value ro 
of the pair distances. The two-body correlation function to find two nucleons at a 
distance r has its maximum around r = 1 f m universally for all light nuclei EKEJ-QIJ 
Thus ro has to be smaller than 1 fm. For the most simple correlated shell model 
ansatz of Eq. (|3-ip or symmetric extensions beyond s-wave and (|3-5[) the exact wave 
function for 3 He and 4 He is to be projected onto the totally antisymmetric spin- 
isospin states X3 an d X4, Eqs. (|2-4p and (I2-5I) . respectively: 



yexact X34 \^ct > _ (6-1) 

For a global adjustment one has to minimize 

3,4 3,4 

/ dv(*i% act - n r(^)? n °^ - r °) ( 6 - 2 ) 

J i=l i<j 

or an extension including higher partial waves but still keep the symmetry in the 
space part. 

The resulting R(u) and R's related to higher partial waves should be independent 
of ro- This requirement should determine the smallest possible value for ro. 
Knowing R(u) one can compare the norms 

3,4 



N exact _ I jy^exactf J| ©(^ _ rQ ) 

J i<j 

3,4 3,4 

N SM = / d y\ HR( Ui )f H0( rij - r ) (6-3) 

i=l i<j 

In addition the short range behavior is not accessible to the single particle picture 
and provides the norm contribution 



jt.Texact 
short 



. 3,4 

/ dVl^rfHOiro-nj) (6-4) 

J i<j 



Finally, one has to keep in mind that only about 90% of the total norm is related to 
the spin-isospin states X3A- The rest is of more complicated structure. Li 3 -* 



§7. Summary 



In nature a nuclear wave function is translationally invariant. Therefore, if a 
shell model picture is a good representation of a nuclear wave function, the single 
particle states have to depend on translationally invariant coordinates. Our choice 
of coordinates Ui = Xi — X relating the individual position vectors x% to the c. m. 
coordinate X fulfills this condition with the additional constraint that they have to 
sum up to zero: Ya=i = 0. Choosing the first n — 1 of them together with the 
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c. m. coordinate one can formulate a shell model Hamiltonian composed of kinetic 
energy containing now also mixed terms V Ui ■ V Uj and single particle potentials 
depending on the coordinates |uj|. Assuming the energetically lowest energy state 
to that Hamiltonian to be a Slater determinant with equal space dependent single 
particle wave functions, R{ui), which is the most simple choice, one obtains nonlinear 
equations for R(ui). They have been worked out for nucleon numbers A = 3 and 4. 
For the special choice of harmonic oscillator potentials the nonlinear equations can 
be analytically solved and that most simple ansatz for the wave function turns out 
to be correct. In the case of general mean field potentials partial wave contributions 
beyond s-states might be necessary. 

In order to shed light on the question how well such a shell model ansatz is 
justified we regarded in some detail three and four nucleons. The corresponding 
shell model Hamiltonian can be written in terms of standard Jacobi coordinates 
and numerically exact solutions can be generated based on the Faddeev-Yakubovsky 
equations. Knowing the exact wave functions one can check the validity of the 
Slater determinant ansatz. Optimization algorithms are provided to perform the 
comparison of exact wave function with the Slater determinant ansatz. 

The main task however, is to confront such a shell model ansatz to realistic three- 
and four- nucleon wave functions (and beyond), which are based on modern two- and 
three-nucleon forces. Clearly at short pair distances the well established repulsive 
nature of the nuclear forces invalidates the shell model ansatz and therefore only for 
pair distances beyond a certain value ro the shell model picture can make sense, if 
at all. To that aim numerical investigations are planned for both, the shell-model 
Hamiltonians and realistic Hamiltonians composed of two- and three-nucleon forces. 

References 

1) B. Giraud, J. C. Hocquenghem, and A. Lumbroso, Proc. Colloque de la Toussuire, lecture 
61 (Fevrier 1971). 

2) D. L. Hill and J. A. Wheeler, Phys. Rev. 89 (1953), 1102. 

3) J. J. Griffin and J. A. Wheeler, Phys. Rev. 108 (1957), 328. 

4) C. W. Wong, Phys. Rep. 15C (1975), 283. 

5) D. Rozpedzik et al., Acta Phys. Polon. B 37 (2006) 2889. 

6) A. Nogga, E. Epelbaum, J. Golak, H. Kamada, H. Witala, D. Rozptjdzik, R. Skibihski, 
and W. Glockle, Four-nucleon force contribution to the binding energy of 4 He, in 19th 
International IUPAP Conference on Few-Body Problems in Physics, Vol. 3 of EPJ Web of 
Conferences, (2010) 05006. 

7) H. Kamada, W. Glockle, Phys. Lett. B 292 (1992), 1. 

8) A. Nogga, H. Kamada, W. Glockle, Nucl. Phys. A 689 (2001), 357c. 

9) A. Nogga, D. Hueber, H. Kamada, W. Glockle, Phys. Lett. B 409 ( 1997), 19. 

10) W. Glockle, H. Kamada, Phys. Rev. Lett. 71 (1993), 971. 

11) A. Nogga, H. Kamada, W. Glockle, Phys. Rev. Lett. 85 (2000), 944. 

12) H. Kamada et al, Phys. Rev. C 64 (2001), 044001. 

13) A. Nogga, H. Kamada, W. Glockle, B. R. Barrett, Phys. Rev. C 65 (2002), 054003. 

14) W. Glockle, H. Kamada, J. Golak, A. Nogga, H. Witala, R. Skibihski, J. Kuros-Zomierczuk, 
Acta Phys. Polon. B 32 (2001) 3053. 

15) W. Glockle, H. Kamada, H. Witala, D. Hiiber, J. Golak, K. Miyagawa, S. Ishikawa, Few 
Body Syst Suppl 8 (1995) 9. 

16) W. Glockle, H. Witala, H. Kamada, D. Hiiber, J. Golak, Few Body Syst Suppl 9 (1995) 
384. 

17) B. Blankleider and R. M. Woloshyn, Phys. Rev. C 29 (1984), 538. 



On a translationally invariant nuclear single particle picture 15 

18) J. L. Friar et al, Phys. Rev. C 42 (1990), 2310. 

19) L. D. Faddeev, Sov. Phys. JETP 12 (1961) 1014. 

20) W. Glockle, The quantum mechanical few - bodyproblem, Springer Verlag 1983. 

21) H. Liu, Ch. Elster, W. Glockle, Computer Physics Communications 147 (2002) 170. 

22) H. Liu, Ch. Elster, W. Glockle, Few-Body Syst. 33 (2003) 241. 

23) A. Stadler, W. Glockle, P. U. Sauer, Phys. Rev. C 44 (1991), 2319. 

24) O. A. Yakubovsky, Sov. J. Nucl. Phys. 5 (1967) 937. 

25) H. Kamada, W. Glockle, Nucl. Phys. A 548 (1992), 205. 

26) M. R. Hadizadeh, S. Bayegan, Few-Body Syst. 40 (2007) 171. 

27) R. B. Wiringa, V. G. J. Stoks, R. Schiavilla, Phys. Rev. C 51 (1995), 38. 

28) R. Machleidt, Phys. Rev. C 63 (2001), 024001. 

29) V. G. J. Stoks et al., Phys. Rev. C 49 (1994), 2950. 

30) E. Epelbaum, Prog. Part. Nucl. Phys 57 (2006) 654 . 

31) S. A. Coon, M. T. Peha, Phys. Rev. C 48 (1993), 2559; J. L. Friar, D. Hiiber, U. van 
Kolck, Phys. Rev. C 59 (1999), 53 ; S. A. Coon, H. K. Han, Few-Body Syst 30 (2001) 131. 



